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ABSTRACT 

In collaboration with NASA/Langley Research Center, a truly 
nonintrusive and nonseeding method for measuring supersonic molecular 
flow parameters had been proposed and developed at Colorado State 
University. The feasibility of this Raman Doppler Velocimetry (RDV) , 
currently operated in a scanning mode, has been demonstrated not only in 
a laboratory environment at Colorado State University but also in a 
major wind tunnel at NASA/Langley Research Center. This report 
summarizes the research progress of this new RDV development. In 
addition, methods of coherent Rayleigh- Brillouin spectroscopy and 
single-pulse coherent Raman spectroscopy are investigated, respectively, 
for measurements of high-pressure and turbulent flows. 
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INTRODUCTION 

Since its first report [Yeh ^nd Cummins, 1964] laser Doppler 
velocimetry (LDV) for measuring fluid flows have under gone more than 
two decades of research, development, and applications; much progress 
has been made. Since the signal strength of LDV depends on scattered 
light intensity from naturally occurring or carefully seeded particulate 
matter in the flow, the accuracy of LDV is limited by the problem of 
particle drag making its effectiveness questionable at supersonic 
speeds. Spectroscopic methods that could probe the motion of atoms or 
molecules directly are therefore of obvious interest. Using laser- 
induced fluorescence with sodium or iodine seeding, speeds of mean 
molecular flow have been determined by measuring either the transit time 
[Proden et al., 1982] or Doppler shift [Miles, 1975; McDaniel et al., 
1983] of the seeding atoms or molecules. The seeding requirement, 
however, is often unacceptable or sometimes difficult to achieve, e.g. 
in combustion or flame. To avoid these difficulties, we have proposed 
the use of coherent Raman spectroscopy for flow velocity measurements 
[She et al., 1981]. The Doppler shift of a fluid component, such as 
nitrogen molecules, is measured directly; with this method, no seeding 
is required whatsoever. 

The feasibility of this ultimate idea, Raman Doppler Velocimetry 
(RDV) has soon been demonstrated in laboratory environments with 
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coherent anti-Stokes spectroscopy (CARS) at Stanford [Gustafson et al., 
1981] and with stimulated Raman gain spectroscopy (SRGS), Inverse Raman 
spectroscopy (IRS), and coherent Stokes Raman spectroscopy (GSRS) by our 
group at Colorado State University [Herring et al., 1981, 1983a; 
HoosmCiller, 1985]. In addition, we have also demonstrated the use of 
IRS technique to measure temperature and pressure [Herring et al . , 

1983a] as well as to map out the two-component velocity distribution of 
a supersonic jet [Moosmuller et al., 1984]. Very recently, using a 
retrometer to minimize the effect of mechanical vibration and to obtain 
simultaneous forward and backward scattering data, the RDV has been 
successfully applied in a major supersonic wind tunnel at NASA/Langley 
Research Center providing measurements of molecular flow velocity, 
temperature, and pressure pertaining to a simple aerodynamic model 
[Exton and Hillard, 1986]. Thus, in six years, the RDV technique is 
becoming an ideal tool for wind tunnel measurements and for aerodynamic 
research at transonic and supersonic speeds. 

In these current versions of RDV, two laser (pump and probe) beams, 
one of which is tunable, are focussed and crossed at the molecular flow 

of interest. As the tunable laser is scanned across the molecular 

resonance, the Doppler shift of Raman resonance, is measured, from 

which the flow velocity v is determined, by 

A^d ” (^p • (1) 

where and are, respectively, the wave vectors of the pump and 
probe laser beams. Since the kinematics of the interaction is quite 
general, the pump and probe method for measuring molecular flow (Eq. 1) 
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should not be restricted to Raman scattering only. Indeed, the use of 
other scattering techniques, such as those for various types of 
conventional laser Doppler velocimetry [She and Kelley, 1985] and a new 
coherent Rayleigh-Brillouin scattering method for molecular flow 
measurement [Herring et al., 1983b] have been established. 

In this report we review the physical principles of molecular light 
scattering and discuss velocity, temperature, and pressure dependences 
in a molecular light scattering spectrum. The methods of extracting 
state parameters of flowing molecular gas from coherent Raman spectra as 
well as other types of coherent light scattering spectra will be 
presented and recent results reviewed. In addition, possible future 
developments as well as the potential and limitation of using coherent 
light scattering spectroscopy for molecular flow measurements will be 
discussed. 

THEORY 

To illustrate the usefulness of light scattering spectroscopy for 
flow measurements, it is perhaps simplest to consider a specific 
example. Figure 1 shows the spectrum of scattered light of simple 
molecules such as nitrogen. The incident laser is assumed to be 
monochromatic. Vibrational and rotational motions of the molecule give 
rise to Raman scattering: vibrational spectrum and associated rotational 
sidebands around 2000 cm and pure rotational spectrum around 
100 cm’^. Near the incident laser frequency, there is a low-frequency 
quasi-elastic scattering spectrum resulting from the translational 
motion of the molecules; its shape depends on the scattering geometry as 
well as the state of the molecular gas. The scattering geometry is 
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determined by the scattering wavevector K which is defined as the 
difference between incident and scattered wavevectors. When the 
scattering wavelength A(2ir/K) is smaller than the mean- free -path L, the 
molecular system exhibits single particle behavior; we have Rayleigh 
scattering with a Doppler broaden spectrum. When A/L » 1, the molecules 
move collectively; scattering from the Doppler shift of this collective 
motion leads to additional peaks at a frequency depending on its (sound) 
velocity; this triplet spectriim is called Rayleigh- Brillouin spectrum. 

Since the location and shape of vibrational Raman, rotational 
Raman, and Rayleigh- Brillouin spectra depend upon pressure, temperature 
and velocity of the flowing molecular gas, each portion of the 
scattering spectrum can be used to determine these flow parameters of 
interest. When molecules are flowing, the peaks of vibrational Raman, 
rotational Raman, and Brillouin spectra, in the non- relativistic limit, 
will be Doppler shifted by an amotint proportional to its flow (drift) 
velocity. Measurement of the Doppler shift provides a determination of 
the flow speed by each spectroscopic technique. Extensive theoretical 
study of Rayleigh-Brillouin scattering spectrum in the past decade 
[Tenti et al., 1974] has made accurate determination of pressure and 
temperature by Rayleigh-Brillouin scattering possible. The random 
collisions of molecules also gives rise to a broadening of rotational 
and vibrational Raman spectra whose widths depend on pressure as well as 
translational temperature. In addition, the Stokes/anti -Stokes 
intensity ratio of vibrational Raman spectra and the relative 
intensities of rotational sidebands can be used to determine, 
respectively, the vibrational and rotational temperatures of the gas. 

For a measurement time greater than the rotational and vibrational life 
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time of the gas, the system Is in thermal equilibrium and the measured 
translational, rotational, and vibrational temperatures should be the 
same. 

Since the spectral resolution of a spectrometer, tjrpically 1 cm’^ - 
30 GHz, is much too coarse for the Doppler shift of interest, <1 GHz, 
the methods of coherent light scattering with resolution limited only by 
laser linewidths must be used. In a coherent light scattering 
experiment, two laser beams, a pump at u>^ and a probe at are crossed 
at the molecules of interest as shown in Fig. 2(a). An exchange of 
energy between the two laser beams via two -photon process takes place 
leading to a gain or loss in the probe beam; in addition, a coherent 
generation (emission) of a third light beam at via four-photon 
process, occurs simultaneously. If the frequency of one laser beam 
either or is tunable, their frequency difference can be made to 
match the characteristic frequency of the meditim, A. Under this 
condition, the scattering process, said to be in resonance, is easily 
detectable giving rise to a peak in the measured spectrum. The 
characteristic (vibration or rotational) frequency of the molecule can 
then be determined this way. There are four tjrpes of coherent 
spectroscopies yielding similar spectroscopic information about the 
medium in question. The choices depend solely on signal-to-noise 
considerations as well as the availability of specific laser systems. 

The schematic descriptions of these processes are shown in Fig. 2(b). 
Both gain and loss spectroscopies are two -photon processes in which the 
probe beam either gains energy (with the probe at the Stokes frequency 
Wp - - A) or loses energy (with the probe at the anti -Stokes 

frequency a> - o>^ - + A) , The emission spectroscopies are the result 
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of four-wave mixing processes generating a third beam at the anti -Stokes 
frequency w - « - + A when the probe is at w or at the Stokes 

6 d X S 

frequency w — w - w- - A when the probe is at w . Although the signal 
e s J> a 

of the two-photon processes is stronger than that of the four-photon 
processes, the signal-to-nolse of these spectroscopies Is In the shot- 
noise limit roughly the same. 

In order to assess the feasibility of a given experiment, a brief 

discussion of signal strength is in order. The signals in question are 

either the change in probe power, AP^, and/or the generated power of the 

third beam, P , respectively for two-photon and/or for four-photon 
a 

processes. The detected signal depends on the laser powers, beam 
crossing angle, and for four-photon processes on the fulfillment of the 
phase -matching condition [Prior, 1983] as well. As the beam crossing 
angle decreases, signal increases at the expense of spatial resolution 
of the measurement. For Gaussian beams, the coherent signals near 
resonance may be calculated from the following formulae: 


lAP 1 - f).G(o) )P 

* p* 2 p p 


for gain or loss spectroscopy (2) 


P^ - >;^(0. 241)G(Wp) Pp for emission spectroscopy (3) 


The exact expression for P^ in terms of gain function with Gaussian 
beams, G(w^), to our knowledge, is not given in the literature. The 
factor 0.241 is derived from a manerical integration under the phase - 
matching condition for four-wave mixing interactions. The gain function 
G, in SI units, may be written as 
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G(Wp) - [4»r^c«^n^/«Wp^]Ng(Wp)(da/dn)P^ (4) 

where N and Pj^ are respectively refractive index at the pump 
frequency, molecular concentration and pump laeor power. For Rnmnn 
processes, gCt^’p) is the Raman lineshape function; it is to be replaced 
by (ftA/kgX)S(K,A) for Rayleigh-Brillouin processes [She et al., 1985]. 
The function S(K,A) is the spectral density function for Rayleigh- 
Brillouin scattering [Tenti et al., 1974], The differential cross 
section do/dO may be obtained from spontaneous scattering measurements 
for both Raman and Rayleigh-Brillouin processes. The signal strength 
can be estimated from Eqs. (2), (3), and (4) for various cases. 

The crossing efficiencies 1^2 *74 two-photon and four- 

photon processes are fvinctions of x - (f/a)tan$ depending on the 
crossing angle 0 , beam radius a, and the focal parameter f. For two 
completely overlapping Gaussian beams, 1 J 2 " ■■ 1. For x larger than 

one, these crossing efficiencies may be calculated approximately from 
the following fonnulae: 


- (2A)tan’^(x^ - 1)*^/^] (5a) 

and 

^2 - (x^ * (5b) 

where x - (f/a)tan^. The crossing efficiencies may be calculated 
exactly by numerical integration by considering beam overlaps. The 
solid curve in Fig. 3 is a plot of the exact for $ - 6.5°. Compared 
to Eq. (5a) shown as the dotted curve, the approximate validity of above 
formulae for x > 1 is evident. 
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The accuracy of any measurement depends on signal -to -noise and 
available measurement time. The accuracy of flow velocity measurement 
by coherent Raman spectroscopy when a laser source is scanned in 
frequency may be expressed in terms of signal- to-noise, measurement time 
and Raman linewidth. A practical formula for estimating measurable 
uncertainty in flow speed is given in the Appendix. 


METHODS OF COHERENT RAMAN SFECTROSCOFT 

After the initial measurement of a subsonic speed in a flow 
using counter -propagating stimulated Raman-gain spectroscopy [Herring 
et al., 1981], supersonic flow in a miniature wind tunnel was made by 
inverse Raman spectroscopy using two overlapping copropagating Gaussian 
beams. The experimental setup is shown in Fig. 4. A cylindrically 
symmetric converging- diverging section (made from Pyrex tubing) with a 
minimum throat diameter of 1,8 mm is followed by a 3. 5 -mm- diameter 
constant-area test section 3 cm in length. A N 2 bottle provides a 
pressure of 550 torr at the input to the wind tunnel while a 
300-liter/min mechanical pump keeps the output section of the test 
region at about 70 torr. This arrangement produces a Mach 2 flow. 

The probe laser is a 5145-A cw single-mode Ar-ion laser. 

Typically, 200 mW of probe laser power is used in the measurement 
region. The high-power ptunp laser is a pulse -amplified cw dye laser; a 
100 mW single-mode cw dye-laser beam is passed through three stages of 
amplification ptjunped by a frequency-doubled NdiYAG laser operating at 
10 pulses/sec. This scheme provides dye-laser pulses of 5-nsec in 
width, 2-MW in peak power, and with a Fourier -limited bandwidth of 
100 MHz. The pump and probe beams are overlapped and directed along the 
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axis of the wind tunnel. Both beams are then focused with a 10*cm lens 
Into the center of the test section. The beam confocal parameter is 
estimated to be 2 mm. After the Interaction, the two beams are 
separated with a prism and filters, and the probe beam Is monitored by a 
high-speed photodiode. To reduce the noise that is due to probe 
intensity fluctuations differential detection on the probe signal is 
employed. The pump dye laser is tuned to around 5845 A, which 
corresponds to Raman resonances from the Q branch of the v - 0 to v - 1 
vibration transition in N 2 . As the dye laser is scanned, the IRS signal 
(after subtraction of the laser noise) is filtered with an RC high-pass 
filter, amplified, averaged by a boxcar integrator, and recorded on an 
x-y recorder. 

To make velocity measurements, the Raman spectra from flowing 
and stationary N 2 compared. From these we measure the Doppler shift 
and hence the velocity. Since the high-velocity and stationary spectra 
are not recorded simultaneously, an absolute -frequency reference is 
necessary to account for possible dye- laser drifts. So a small part of 
the dye laser beam is sent into and I 2 absorption cell, and the Doppler- 
free saturated- absorption spectrum of I 2 is recorded simultaneously with 
each IRS spectrum. This permits comparison of Raman spectra obtained at 
different times by simply overlapping the I 2 spectra. The many lines 
recorded in an I 2 spectrum correspond to hyperfine components of a 
single rotational -vibrational transition near 5854 A. The results of a 
typical 10-min velocity measurement in the copropagating case are shown 
in Fig. 5. The pump and probe powers in the interaction region are 
1.5 MW and 200 mW, respectively. The shift of 110 MHz in the Raman 
spectrvim corresponds to a flow velocity of 475 m/sec. The precision of 
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our measurements is better than 5%. Twelve different measurements were 
made for the same flow conditions. This group gives an averaged Doppler 
shift of 109 MHz with a standard deviation of 3 MHz, corresponding to a 
flow velocity of 471±13 m/sec. 

For the same flow, a sample temperature measurement is shown in 
Fig. 6, which shows two differept 30-GHz scans over the lines J - 8, 9, 
10. The top curve illustrates the relative intensities of these lines 
at room temperature, and the bottom curve shows the intensities for 
supersonic flow with a clear reduction in the J - 10 line, indicating 
expansion cooling. Since the linewidths for different J values in a 
given flow condition are assumed to be the same, the ratio of the 
heights of the J - 8 and J » 10 lines can be used to determine the gas 
temperature. The bottom trace gives a temperature of about 170 K. For 
three consecutive measurements, the maximxim variation in the measured 
temperature is about 30 K for room- temperature measurements but is only 
10 K for the 175 -K measurements. The smaller absolute errors at lower 
temperatures are expected for a given signal -to -noise ratio. 

Density measurements can also be made from the traces in Fig. 6, 

The density at a room- temperature , 70-torr, stationary condition can be 

determined from the ideal gas law. With this value, the ratio of the 

integrated Raman intensities is used to determine the density under 

the supersonic flow condition. Using data (J - 8) shown in Fig. 6, we 

18 - 3 

obtain a density of about 2.2x10 cm . As with the temperature 
measurements, the results of three consecutive measurements on the same 
flow Indicate uncertainties of about 10% or less. For density 
measurements we compare integrated intensities rather than peak heights 
because the linewidths (for the same J value) for different flow 


11 


conditions ar« expected to have different Doppler and pressure 
broadening. In fact, the flow velocities in our Interaction region may 
not be uniform, which could lead to a small but detectable additional 
broadening. In principle, if the temperature, laser linewidth, and this 
additional (velocity) broadening are known, the Raman llnewidths can be 
used to measure the pressure of the supersonic flow directly. 

Using a retrometer system and overlapping Gaussian beams, Exton and 
Hillard were able to obtain both copropagating (forward) and counter- 
propagating (backward) inverse Raman spectra of the molecular flow In 
the same frequency scan. Such a combined spectrum is schematically 
shown in Fig. 7. The difference between forward and backward Doppler 
shifts, W, may be related to the velocity component along the beams, 

W - (dWjj)j - “ 2Wp(v/c)cos(jp (6) 

where v is the magnitude of the velocity vector, the angle between the 
(colinear) laser beams and the velocity vector, and u>^ is the probe 
frequency. Since the pressure broadening contributes to linewidth 
equally in both forward and backward geometi^, while Doppler broading 
(temperature dependent) contributes much less in the forward scattering, 
the measurement of both forward and backward llnewidths allows a 
determination of both pressure and temperature of the flowing gas. A 
chart showing the temperature and pressure dependencies of forward and 
backward llnewidths is shown in Fig. 8. Using this unified approach, 
Exton and Hillard have been able to measure distributions of molecular 
flow velocity, static pressure and translational temperature in a 
supersonic wind tunnel in both free-stream operation and near a simple 
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aerodynamic model [Exton and Hillard, 1986] demonstrating the 
feasibility of RDV for wind tunnel applications. 

If two probe beams with proper crossing angles to the pump beam are 
used, two velocity components of a flow can be measured simultaneously 
[HoosmClller et al., 1984], To determine two perpendicular velocity 
components we have to make and ) perpendicular to 

each other. For a given Raman shift A and probe -laser wavelength A , 
this leads to a unique geometry for the wave vectors, as shown in Fig. 9 
for A — 2330 cm ^ and A^ - 514.5 nm. As can be seen in this figure, it 
is possible to measure directly two perpendicular velocity components by 
crossing the two probe beams with an angle of only 6.5” to the pump 
beam. The spatial resolution is determined by the minimum spot size and 
the crossing angle. In our case, a spot size of 35 pm and a crossing 
angle of 6.5” results in an interaction length of 0.7 mm, which 
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corresponds to an interaction volume of 0.003 mm . The smaller 
interaction length also results in a smaller signal compared with the 
colinear arrangements. In our case the signal is reduced by roughly a 
factor of 30. 

The experimental set up is shown in Fig. 10. The supersonic jet is 
generated by expanding N 2 from a back pressure of 9.8 bars through a 
nozzle with a diameter D of 1.51 mm into a background pressure of 
53 mbars. A N 2 bottle provides gas with sufficient pressure, and the 
background pressure is maintained by a 4300-llter/min mechanical pump. 
Using an empirical formula [Ashkenas and Sherman, 1966], we determine 
the length of the supersonic region to be 13.7 mm. The pump and probe 
beams are as described before and they are focused by lenses with 20-cm 
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focal length. At their foci the probe beams intersect symmetrically the 
focus of the pump beam at an angle of 6.5”. 

To make velocity measurements, Raman spectra of stationary and 
flowing N 2 have to be compared. Initially we attempted to record the 
spectrum of stationary in a copropagatlng-beam arrangement 
simultaneously with the spectrum of flowing which is taken in a 
crossed-beam arrangement. Comparing spectra of stationary N 2 with both 
arrangements, we observed a systematic frequency shift on the order of 
100 MHz. This shift is possibly due to a Stark effect and will be 
investigated further in the future. To avoid this complication, we have 
taken reference spectra of stationary N 2 with the setup used for the 
flow measurement. To account for a possible frequency drift of the dye 
laser between the different scans, we use part of the cw dye -laser beam 
to record a saturated absorption spectrum simultaneously in I 2 . The 
hyperfine components of an I 2 line are then used as an absolute- 
frequency reference. The result of a typical 10-min measurement of two 
velocity components on the axis of the jet is shown in Fig. 11. Both 
peaks are shifted nearly symmetrically relative to the peak of a 
spectrum of stationary N 2 . This Indicates a flow velocity perpendicular 
to the pump beam, as expected. We have mapped out the velocity 
distribution of a center plane of the jet sequentially. The results of 
these measurements are shown in Fig. 12, along with the boundary of the 
supersonic jet determined from a simple Schlleren setup. The velocity 
vectors behave in general as expected. We observe expansion and 
acceleration with increasing distance from the nozzle exit. We have 
also Investigated the temperature distribution within the Jet. 
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Because of the need to frequency scan one laser beam to obtain the 
desired Raman spectrvim for determining the Doppler shift, a velocity 
measurement as described above typically requires between 5 to 10 min. 
Althou^ the relatively long measurement time of the current RDV Is 
acceptable In certain wind tunnel applications, it Is a severe 
limitation for tunnels with short run-time or tunnels under turbulent 
flow conditions. This fact has been recognized even in our initial 
proposal [She et al., 1981] which contains some suggestions to reduce 
the measurement time. More recently, the possibility for velocity 
measurement with single-pulse coherent Raman spectroscopy has been re- 
examined [She, 1983]; It concluded that with the development of a 
specialized probe laser and high- resolution spectrum analyzer, single- 
pulse flow velocity measurements should be possible. Very recently we 
have conceived yet another method for performing single -pulse velocity 
measurement with coherent Raman spectroscopy [She, 1987]. We suggested 
the use of a narrow-band pulse pump laser together with a cw probe laser 
consisting of several discrete frequency components separated by an 
amount smaller than the Raman linewidth of Interest. In this case, a 
discrete spectrum of coherent line emissions (CARS or CSRS) , one 
associated with each probe frequency component, can be detected. 

Analysis of this discrete spectrum by a suitable interferometer should 
allow the determination of the spectral peak, Doppler shift and hence 
the flow velocity. The advantage of this new technique over the other 
single-pulse methods lie in the fact that the required probe laser can 
be implemented from a cw single -frequency laser by the method of 
electro-optic (either amplitude or frequency) modulation. With the 
successful implementation of any single-pulse measurement method, the 
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RDV technique should become a truly noninvaslve, real-time probe for 
aerod 3 mamlc research under any flow conditions imaginable, laminar or 
turbulent . 

As described in the Appendix, the accuracy of velocity 

determination is limited by the available signal -to-noise in an 

experiment. With current scanning method, the nominal speed uncertainty 

is around 1 m/sec; the uncertainty Is expected to increase for the 

proposed single-pulse measurement for lack of averaging. In principle, 

with higher pump power the coherent Raman signal can be made to 

Increase. However, the presence of optical Stark effect, as pointed out 

previously, will be a fundamental limit that cannot be ignored. Since 

-1 2 

Stark shift which is about 0.1 cm per 1 TW/cm [Rahn et. al., 1980] 
depends only on intensity, higher pump power may be used in a life-size 
wind tunnel without over penalizing spatial resolution. Optical Stark 
effect is of interest in its own right in nonlinear optics; its 
continued study is in progress. 

COHERENT RAYUSIGH-BRILLOUIN TECHNiq^ 

As mentioned in the Introduction, both Raman and fluorescence 
techniques are good spectroscopic methods for measuring molecular flows; 
however, they have a common problem. The pressure broadening of the 
Q-branch (v — 0 -► v — 1) Raman transition in molecular nitrogen is 
2.4 GHz/atm. This broadened linewidth makes the determination of the 
Doppler shift more difficult and the technique less sensitive for 
velocity measurements at higher pressures. The pressure effect on 
fluorescence is even more destructive. In addition to line broadening, 
the fluorescence of iodine and sodium is quenched by one order of 



magnitude at a nitrogen pressure of 0.1 atm. Thus these spectroscopic 
techniques for velocity measurements are becoming Increasingly 
Impractical at pressures exceeding 1 atm. Although most wind tunnels 
use pressures lower than 1 atm, there are important large facilities 
that are operated at higher pressures. An outstanding example is the 
new NASA National Transonic Facility, which typically operates at a flow 
from Mach 0.2 to 1.2 and a nitrogen pressure between 1 and 9 atm. For 
velocity measurements of such high-pressure flows, a different 
spectroscopic technique is apparently needed. We now discuss the use of 
stimulated Rayleigh-Brillouin-gain spectroscopy for measuring atomic- or 
molecular -flow velocities. This technique is most suitable for flows 
with a pressure greater than 1 atm because, uni ike the Raman or 
fluorescence techniques, its sensitivity improves as the gas pressure 
increases . 

Rayleigh-Brillouin-gain spectroscopy is a form of coherent 
scattering spectroscopy in which the frequencies of the probe, and 
pump, lasers, one of which is tunable, are within a few GHz of each 

other. Similarly to stimulated Raman scattering, the Rayleigh-Brillouin 
spectrum shows gain on the stokes side > u>^) and loss on the anti- 
Stokes side (o >2 < , leasing to a zero crossing at " ^2* Figure 13 

shows the qualitative difference of stimulated Rayleigh-Brillouin-gain 
spectra of nitrogen at 0.9 and 4.1 atm. The fact that the signal-to- 
noise ratio of the 4.1-atm spectrum is much better and its width 
narrower is evident. Since the peak gain increases quadratically at 
high pressure because of the narrowing of Brillouin linewidth as the 
pressure increases, the method of stimulated Rayleigh-Brillouin-gain 
spectroscopy works better for higher pressures. 
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Similarly to the Raman technique » the flow velocity v may be 
determined from the measured Doppler shift 

where and are, respectively, the wave vectors of the pxjunp and 
probe lasers. In addition to the flow velocity, fitting the Rayleigh- 
Brillouin-gain spectrum to a theoretical model can yield a determination 
of the local values of thermal conductivity, viscosity, temperature, and 
gas density. 

The experimental setup is shown in Fig. 14. The probe laser is a 

single-mode cw dye laser (l-MHz linewidth) , and the pump laser consists 

of an identical cw dye laser along with a three-stage pulsed amplifier 

(peak power 1 MW, 5-nsec pulses, 10 pulses/sec, and 150-MHz linewidth). 

The two counterpropagating laser beams were crossed at an angle of 1® 

and focused with 20-cm lenses to a waist of approximately 25 and an 

interaction length of 3 mm. The flow from a converging nozzle was 

directed at an angle of 15® with respect to the crossed laser beams. 

This measurement was made in the turbulent subsonic region (after the 

normal shock) about 3 cm downstream from the nozzle. The 1-mm-diameter 

3 

nozzle was backed with 5.69x10 torr of nitrogen and exhausted into room 
air. After it had traversed the interaction region, the probe beam was 
monitored with a high-speed photodiode followed by a high-pass filter 
and an amplifier. The signal was then averaged by a boxcar integrator. 
As the pump laser is frequency scanned across the probe -laser frequency, 
a stimulated Rayleigh-Brillouin-gain spectrum is recorded. The gain 
spectra with flows directed in opposite directions and with stationary 
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air are compared, and the Doppler shifts that are due to the flow are 
determined. The result for flow- velocity measurements using stimulated 
Raylelgh-Brillouin-gain spectroscopy is shown in Fig. 15, where spectra 
of flow directed along the pump and probe beams are shown. The spectrum 
of stationary room air (not shown), which has a comparable signal-to- 
nolse ratio, would fall between the Rayleigh-Brillouin spectra, if 
superimposed. An iodine fluorescence spectrum was recorded 
simultaneously with each Rayleigh-Brillouin-gain spectrum for a 
frequency reference. The shift between peaks (or zero crossing) is 
determined to be 1.2 GHz, corresponding to a relative velocity of 
350 m/sec. The gain spectra were repeated three times for each 
direction, yielding nine different measurements of the same flow. These 
measurements had a maximxjkm discrepancy of 60 m/sec. Since this 
discrepancy is Independent of the flow velocity, the percentage 
measurement uncertainty would be smaller for faster flows. 

The technique is obviously applicable to supersonic flows. 

Although the measurements reported here are done at atmospheric 
pressure, the Rayleigh-Brillouin signal should increase drastically with 
gas pressure. Flow-velocity measurements and the determination of other 
local flow parameters (viscosity, temperature, and density) at pressures 
higher than 1 atm can be made nonintrusively with lasers using 
stimulated Rayleigh-Brillouin-gain spectroscopy. Other advantages in 
Rayleigh-Brillouin spectroscopy for flow measurements include the facts 
that the same laser system may be used for any molecular system and 
that, unlike vibrational Raman scattering, it works for atomic systems 
as well. Like coherent Raman spectroscopy, other forms of coherent 
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Raylel^-Brlllouln spectroscopies at different crossing angles can also 
be used to probe atomic and molecular flows. 

CONCUJSI(» 

In this report, we have presented a unified discussion of coherent 
light scattering spectroscopy and its application for a nonintruslve 
determination of the state of flowing molecular gas. These include 
measurements of flow velocity, gas temperature and static pressure. The 
feasibility of using coherent Raman spectroscopy for this purpose has 
been demonstrated not only in the laboratory but also in a major wind 
tunnel. The advantage of using coherent Rayleigh- Brlllouln scattering 
for measuring high pressure flow Id discussed and Its feasibility 
demonstrated In the laboratory. In addition, the possibility and 
methods of single-pulse flow measurements Is discussed; Its successful 
Implementation will make the RDV method a truly nonintruslve real-time 
probe In aerodynamic research for measuring any flow Imaginable, laminar 
or turbulent. 
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Fig. 1. Light scattering spectrum of diatomic molecule, e.g. nitrogen. 


Coherent Raman processes: (a) pump and probe beams crossing at 

an angle $ . (b) energy diagram describing four types of 

processes . 
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Fig. 3. Crossing efficiency of the two-photon process as a function of 
f/a (6 - 6.5*). 
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Fig. 4. Experimental setup for inverse Raman spectroscopy using 

copropagating geometry. [Following Herring et al., 1983a.] 
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Sample copropagating frequency shift of J - 8 for a 475 m/sec 
flow. [Following Herring et al., 1983a.] 
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Sample temperature measurements showing the relative 
intensities of J - 10 (right-most peak). [Following Herring 
et al. , 1983a. ] 
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Sketch of the shift and breadth of the forward and backward 
scattered Raman lines relative to the unshifted (zero velocity) 
position. [Following Exton and Hillard, 1986.] 
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Backward and forward Voight line breadths (FWHM) are shown for 
the N 2 Q-branch (J - 10) line over a wide range of temperature 
and pressure. [Following Exton and Hillard, 1986.] 
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Fig. 9. Geomet^ of the wave vectors for pump beam (L) and probe beams 
I (k . k ). [Following Moosmaller et al., 19847] 
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Fig. 10. Experimental setup for the measurement of two perpendicular 
velocity components with IRS. [Following Moosmdller et al., 
1984.) 
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shifts of the two peaks correspond to two velocity components 
[Following Moosiiller et al., 1984.] 



Fig. 12. Flow- velocity and temperature distribution in the supersonic 
jet. [Following Moosmuller et al., 1984.] 
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Fig. 13. Typical stimulated Rayleigh-Brillouin-gain spectra of nitrogen 
at 0.9 and 4.1 atm, [Following Herring et al., 1983b.] 



Fig. 14. Experimental setup for flow velocity measurements with 

stimulated Rayleigh-Brillouin-gain spectroscopy. [Following 
Herring et al., 1983b.] 
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Fig. 15. Typical result of flow- velocity measurements. The obvious 
shift between the peaks of these Rayleigh-Brillouin-gain 
spectra corresponds to a relative velocity of 350 m/sec. 
[Following Herring et al., 1983b.] 


! 


I 

i 

j 




31 


APPENDIX 

Thicertalnty in Spectral Peak Deterainatlon 

In looking for a new spectroscopic feature or carrying out a 
precision measurement with tunable lasers, one Is concerned with the 
availability of adeqviate signal -to -noise for the experiment. Once 
proceeding with the experiment, one simply uses as much measurement time 
as needed to achieve the best possible relative accuracy and, if 
necessary, by repeating the measurement many times. The absolute 
accuracy of the experiment is usually limited by the availability of a 
standard reference, nearby in frequency, for comparison. In applied 
spectroscopy, one usually deals with a known feature whose spectrum has 
been well characterized under a given environment, e.g. temperature, 
pressure, etc. Probing a sli^t change or shift of the spectroscopic 
feature is used as a tool to determine the experimental condition in 
question. Because a dynamic change is of interest the time, with which 
a spectroscopic measurement is made, is thus limited. In this case the 
relative measurement accuracy depends on available measurement time as 
well as signal -to -noise ratio (S/N) . Since the changes are being 
measured, the absolute accuracy in frequency is usually not an issue. 

The dependence of measurement accuracy on S/N and measurement time 
is derived and discussed in this Appendix. Although only coherent Raman 
spectroscopy and the associated achievable accuracy for the flow 
velocity measurement will be considered, the analysis should have 
general validity for any type of spectroscopy when a laser source is 
scanned in frequency. Since a pulsed laser is used, the signal is 
detected at a repetition rate r with an integration time constant t by a 
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box car Integrator. Each detection Interval At, which is set by the 
pulse width or box car window, is determined by the transform- limited 
bandwidth of the detection system. Along with these quantities, the 
signal- to-nolse in a detection Interval, S/N, which depends on laser 
power and its Intensity fluctuation, d, and the measurement time, T, are 
the only parameters input to the analysis. The aim of the analysis Is 
to provide a working guide for estimating the achievable accuracy of an 
anticipated measurement, not to give a statistical treatment of the 
measurement. Thus, the result of the analysis is useful but not 
absolute; it is subjected to the interpretation of detection bandwidth, 
1/At, and integration time constant, t. 

For coherent Raman spectroscopy, the slgnal-to-nolse in a detection 
Interval is well known. It is given, respectively for SRGS (IRS) and 
CARS (CSRS) as 


S/N - 


[8^ + 2ho> /QP At]^/^ 


and S/N - 


+ 2ttw /QqP At]^/^ 


where G, A, <a , o> , Q, P, and r/ are, respectively, Raman gain, Planck 
constant, Stokes frequency, anti-stokes frequency, detector quantum 
efficiency, probe power, and a reduction factor of the CARS (CSRS) 
process which is due to mismatch of optical phases and beam profiles. 

The laser intensity fluctuation, 5, is the fractional uncertainty of the 
probe laser for the SRGS (IRS) process and of twice of the pump plus 
that of the probe laser for the CARS (CSRS) process. It is worthwhile 
to note that in the shot-noise limit SRGS and CARS have comparable S/N. 
However, if the experiment is limited by laser intensity fluctuation 
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noise, CARS has more favorable S/N xmless the stimulated Raman gain, G, 
is nearly unity. In an integration time constant t, there are n (-rt) 
detection intervals averaged to give a slgnal-to-noise of the 
measurement , 


S/H - (S/N)7S - (S/N)7(it)‘ (Al) 

For simplicity, a symmetric bell-shape spectrum is assumed. The 
uncertainty in determining the spectral peak obviously depends on S/W of 
the measured spectrum. More Specifically, the following Gaussian and 
Lorenzian spectra are considered: 

y - yo«*Pt*(»' - »^q)V 2^^] and y - Yq/Ki^ - »/q)^ + 0^] 

For both these spectra, the maximum spectral change occurs at 
t/ - i'q “ ±/3, where the spectrxun is reduced from the peak by 0.50 and 
0.61, respectively, for Gaussian and Lorenzian profiles. Here, the 
uncertainty in measurable frequency, Aw, may be related to the 
fractional spectral uncertainty as 

dy/du - y//3 or Ai/ - (dy/y)/3 - (A2) 

One may locate the spectral peak by locating the points with maximum 
slope and then declare their midpoint as the peak. For all intents and 
purposes, the points with maximum slope may be regarded as the half- 
power points. The S/W in Eq. (Al) refers to the signal-to-noise at 
these frequencies. The optimum situation is to scan the laser at a 
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rat«, R, which covers the minimum resolvable spectrum Av In an 
Integration time constant, t, l.e. R - Ai//t. The adeqtjuite measurement 
time Is somewhat arbitrary, although It Is quite reasonable to set 
T - 4i8/R. These conditions together with Eqs. (Al) and (A2) yield the 
measurement uncertainty and Integration time constant; they may be 
expressed in terms of repetition rate r, measurement time T, and signal- 
to>nolse S/N as 

Av - (4/rT)^/^/9(S/N)"^/^ and t - (A3) 

Eqxiatlon (A3) is the main result of the present analysis. For example, 

consider an experiment with the gain at the half -power point G - 0.05, 

T - 300 sec, p - 250 MHz, r - 10, Q - 0.1, w - 5.10^^ Hz, 

s 

-8 

At - 10 sec, and — 50 mW in the shot-noise limit. For the peak 
rotational line (J 8) In the Q-branch nitrogen Raman spectrum, this 

requires a pump power less than 0.5 MH. The calculated values are 
S/N - 434, A« 0.48 MHz, and t - 0.14 sec. If the experiment is 
limited by laser intensity noise, S/N is much smaller, e.g. 10. In this 
case, one calculates Ai/ - 6.13 MHz and t — 1.78 sec. The measured 
signal- to-noise of the spectrum for these cases should be, respectively, 
S/U - 514 and 42. The latter value has been experimentally achieved 
routinely. 

For a given accuracy in Ai/, the measurable uncertainty in flow 
speed depends on the geometry of the experiment. If the Raman line is 
pressure broadened, the llnewidth /3 remains constant and independent of 
the experimental geometry. The accuracy in speed measurement would be 
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Av — 2wA»//Ak, where Ak - |ic. - Ic | (A4) 

X s 

Continuing on the above examples, the extreme values for Ak are 
approximately 1.5*10^ m’^ (copropagating geometry) and 2.0»10^ 

(counter -propagating geometry). In the shot-noise limit, these 
correspond to a speed accuracy of 2 m/sec and 0.15 m/sec. For the 
laser-noise limited case, the measurable accuracies in velocity are 
26 m/sec and 2 m/sec. If the spectrum is partially Doppler broadened, 
the linewidth should be wider for the counter -propagating geometry, and 
the uncertainty in velocity should be correspondingly worse than the 
values 0.15 m/sec and 2 m/sec given above for the counter -propagating 
geometry. In any case, with the pump power less than 1 MW, pump and 
probe beams overlapping, one may expect a typical velocity uncertainty 
within 1 m/sec in the shot-noise limit and 15 m/sec with laser intensity 
noise , 

In the interest of spatial resolution, the probe and pump beams may 
be made to cross at an angle. This may decrease the Raman gain by a 
factor up to 10 or even 100. Because of the -2/3 power dependence of Ai/ 
on S/N, the uncertainty in measurable frequency and velocity increases 
only by a factor up to 5 or 22 correspondingly. Even in the laser-noise 
limited case, these frequency uncertainties are still much less than the 
Raman linewidth, suggesting the validity of the measurement technique in 
the crossed geometry as well. Furthermore, this analysis remains valid 
for the three -component velocity measurements provided three times probe 
power (150 mW for the above example) is used. Although a pulsed pump 
laser with a repetition rate r was considered, Eq. (A3) is also valid 
for continuous -wave scanning laser spectroscopy if r is replaced by 
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1/At. In this case, phase 'sensitive detection method with the modulated 
signal contained in the frequency band 1/At, should be used. 


